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I. INTRODUCTION 


Precise determinations of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements 1] 
are important to check the consistency of the Standard Model and search for new physics. 
The value of sin(20!), where 0! is one of the angles of the Unitarity Triangle 2] is now 
measured with high precision: sin(2</> 1 ) = 0.731 ± 0.056 |3|. However, this measurement 
contains an intrinsic ambiguity: 20! <—> 7 r —20x. Various methods to resolve this ambiguity 
have been introduced p[, but they require very large amounts of data (some impressive first 
results notwithstanding p]). 

We suggest a new technique based on the analysis of B° —> Dh°, followed by the three- 
body decay of the neutral D meson, D —> Ks^ + ^~ ■ Here we use h° to denote a light 
neutral meson, such as 7T°, r], p°, u. The modes B° —> Dcph 0 , utilizing the same B decay 
but requiring the D meson to be reconstructed via CP eigenstates, have previously been 
proposed as “gold-plated” modes to search for new physics effects j:6j. Such effects may 
result in deviations from the Standard Model prediction that CP violation effects in b —> 
cud transitions should be very similar to those observed in b —> ccs transitions, such as 
B° —> J/'ipKs- Detailed considerations have shown that the contributions from b —> ucd 
amplitudes, which are suppressed by a factor of approximately 0.02 '7|, can be taken into 
account [§]. Consequently, within the Standard Model, studies of B° —> D C ph° can give 
a measurement of sin(20!) that is more theoretically clean than that from B° —» J/^Kg- 
However, these measurements still suffer from the ambiguity mentioned above. 

In the case that the neutral D meson produced in B° —> Dh° is reconstructed in a 
multibody decay mode, with known decay model, the interference between the contribut¬ 
ing amplitudes allows direct sensitivity to the phases. Thus 20i, rather than sin(20i) is 

extracted, and the ambiguity 20i <-> n — 20i can be resolved. This method is similar 

to that used to extract 03, using B ± —> DK ± followed by multibody D decay f.9, |l0j; in 
the 03 analysis the ambiguities in the result are also reduced compared to more traditional 
techniques ]l ll. 0 . In addition, with current B factory statistics, better precision on 0 3 is 
obtained using multibody D decay. 

There are a large number of different final states to which this method can be applied. 
In addition to the possibilities for h°, and the various different multibody D decays which 
can be used, the method can also be applied to B° —> D*h°. In this case, the usual care 
must be taken to distinguish between the decays D* —> Dir 0 and D* —> D'y 0 . 

In this paper we concentrate primarily on the decay B° —> Dn° with D —> (and 

denote the decay chain as B° —> (Ks7t + ti~) D vr°). This multibody D decay has been shown, 
in the 03 analysis, to be particularly suitable for Dalitz plot studies. In the remainder of 
the paper, we first give an overview of the relevant formalism, and then turn our attention 
to Monte Carlo simulation studies of B° —> (K S TT + 7T~) D h°. We attempt to include all 
experimental effects, such as background, resolution, flavour tagging, and so on, in order to 
test the feasibility of the method. Based on these studies, we estimate the precision with 
which 20i can be extracted with the current B factory statistics. 

Before turning to the details, we note that this method can also be applied to other 
neutral B meson decays with a neutral D meson in the final state. In particular, the 
decay B° —>■ DKg has contributions from b —> cus and b —> ucs amplitudes, which have a 
relative weak phase difference of 0 3 . Thus a time-dependent Dalitz plot analysis of B° —> 
(A 57 r + 7 r _ ) D K s can be used to simultaneously measure 0! and 0 3 |l4(, and can test the 
Standard Model prediction that CP violation effects in b —> cus transitions should be, to a 
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good approximation, the same as those in b —> ccs transitions. Furthermore, modes such as 
Bg —> (/i5-7r + 7r _ ) D 0 can in principle be used to measure the weak phase in B®-B° s mixing. 
However, our feasibility study is not relevant to B ° decay modes, which cannot be studied 
at a B factory operating at the T(4S') resonance. 


II. DESCRIPTION OF THE METHOD 


Consider a neutral B meson, which is known to be B° at time t t ag (for experiments 
operating at the T(4S') resonance, such knowledge is provided by tagging the flavour of the 


other B meson in the T(4Sb 
the B meson is given by [l5j 


BB event). At another time t s i g the amplitude content of 


£°(At)) = e - |At|/2T s° (Jh°) cos(AmAf/2) - i- |B°) sin(AmAt/2) j , (1) 

where At = t S [ g —t tag , r^o is the average lifetime of the B° meson, Am, p and q are parameters 
of B° — B° mixing (Am gives the frequency of B° — B° oscillations, while the eigenstates of 
the effective Hamiltonian in the B° — B° system are | B±) = p \B°) ± q B °^), and we have 

assumed CPT invariance and neglected terms related to the B° — B° lifetime difference 0. 
In the following we drop the terms of e I a *I/ 2t bo. 



FIG. 1: 


Diagram for the dominant colour-suppressed amplitude for B° —> Dtt°. 


Let us now consider the decays of the B meson to Dh°. At first, we consider only the 
favoured b —> cud (and charge conjugate) amplitude, shown in Fig. U1 Then the D meson 
produced by B decay at time t sig is given by the following admixture: 


D b° (At)) — 


ZF 


cos(AmAt/2) — i-r] h o(—l) 1 
Q 


D°) sin(AmAt/2), 


( 2 ) 


where we use r^o to denote the CP eigenvalue of h°, and l gives the orbital angular momen¬ 
tum in the Dh° system [17J . 

The next step is the multibody decay of the D meson. We use D —> Ks7r + n~ for 
illustration. We follow jp] and describe the amplitude for a D° decay to this final state as 
/(m^_, m 2 ), where m 2 + and rn 2 _ are the squares of two body invariant masses of the Kg^ + 
and Ks^~ combinations. Assuming no CP violation in the neutral D meson system, the 
amplitude for a D° decay is then given by /(m 2 ,m^_). The amplitude for the B decay at 
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time t S i g is then given by 

Mgo(Ai) = /(m^, m+) cos(AmAf/2) — i-rj h o(— l)*/( m +, m 2 _) sin(AmAf/2). (3) 

Similar expressions for a state which is known to be B° at time i tag are obtained by inter¬ 
changing B° * -> B° , D° < > D°, p < -> q and m 2 + <-> m 2 _\ 


B°(At)) 

= e 

-\At\/2 r B o 

Deo (At)) 

-( 

D cos 

M b o (At) 

= f(m 2 + ,m 2 _ 


B 0 ^ cos(AmAt/2) — i— B 0 ^ sin(AmAf/2) j , 


D°) cos(AmAf/2) — i-rj h o(— l) z sin(AmAf/2) 
s(AmAZ/2) — i-r} h o(—l) l f(m 2 _, m 2 + ) sin 

p 


)■ 


( 4 ) 

( 5 ) 

( 6 ) 


In the Standard Model, |g/p| = 1 to a good approximation, and, in the usual phase 
convention, arg (q/p) = 2<pi. Then 

M-jg 0 (At) = /(mi, m+) cos(AmAZ/2) — ie~ l2(t>1 r] h o(— l) z /(m+, m 2 _) sin(AmAf/2), (7) 
M B o(At) = f(m+,m 2 _)cos(AmAt/2)—ie + ' l2<l>1 ri h o(—l) l f(m 2 _,m+)sm(AmAt/2), (8) 

and it can be seen that once the model f (m 2 + ,m 2 _) is fixed, the phase 2<f>i can be extracted 
from a time-dependent Dalitz plot fit to B° and B° data. 

At this point it is instructive to compare to the B ± —»• DK ± analysis jg]. In that case 
we obtained time-independent expressions 

M b - = f(m 2 _,m 2 + )+r DK e^ SDK ^ 3) f(m 2 + ,m 2 _), (9) 

M b + = f(m 2 +1 m 2 _) + r DK e l< ' 5DK+rf ’ 3 ' 1 f(m 2 _ : m 2 + ), (10) 


where Tdk is the ratio of the magnitudes of the contributing (suppressed and favoured) 
decay amplitudes (tdk = A(B~ —> D°K~) /A(B~ —> D°K~ ) ), and 5dk is the strong phase 
between them. It can be seen that the role of r ok in the time-independent analysis is 
taken by the expression tan(AmAf/2) in the time-dependent case. Furthermore, in the 
time-dependent case, there is no non-trivial strong phase difference. Therefore, the time- 
dependent analysis has the advantage that there is only one unknown parameter, which 
partly compensates for the experimental disadvantages that are accrued. 

We now consider the effect of the Cabibbo-suppressed b —> ucd amplitudes, shown in 
Fig. El The magnitude of this amplitude is expected to be smaller than the Cabibbo-favoured 
diagram (Fig. QJ by a factor of 


TDtt° 


A(B° - 



v ub v; d 

A(B° - 

-> D°7T°) 


v cb v: d 


( 11 ) 


Since this simple approximation neglects hadronic factors, it is the same for all h°, though the 
precise values will depend on the final state. We denote the strong phase difference between 
the two amplitudes as 8o n o (which, in general, will be different for each h °). Including this 
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FIG. 2: Diagram for the colour- and Cabibbo-suppressed amplitude for B a —> Dir 0 . 


amplitude, the expressions Eqs. d and |H] are replaced by 

Mgo(Ai) = cos(AmAf/2) 

ie~ l2< t ,1 ‘q hQ (—\) 1 f(m 2 + ,m 2 _) + r D ^oe l ^ D ^ 0+ ^ f(m 2 + ,m 2 _) 

M B o(At) = f{m 2 + , m 2 _) + r Dn oe v ^DnO+fo) f(m 2 + , m 2 _) cos(AmAt/2) — 

i e +‘ l2< t > i^ 0 (_l) 2 f{m 2 _,m 2 + ) + r Dn ae^ 5D ^°-^ f(m 2 + , m 2 _) sin(AmAf/2). 


( 12 ) 

(13) 


In principle, therefore, it is possible to extract all four unknown parameters 
(20!, r Dn o, 5 dtt°, 03) from the time-dependence of the Dalitz plot. However, due to the small¬ 
ness of r Bn o, this is highly impractical. On the other hand, the above formulation allows 
us to generate simulated data including the suppressed contribution, and thus estimate the 
effect of its neglect. 

Again, we note that in the case of B° —> DK$, the ratio of amplitudes is not small 
C r DK s 0.4), and in this case both 0i and 03 can be extracted from a time-dependent Dalitz 
plot analysis. In fact, the size of r B K s nrakes this mode quite attractive for the measurement 
of 03. Similarly, a time-independent analysis can be performed using B° —> DK *°, with 
K*° —► K “7r + , but in this case additional uncertainty arises due to possible contributions 
from nonresonant B° —> DK~tt + 0- 


III. FEASIBILITY STUDY 

The potential accuracy of the 0! determination is estimated using a Monte Carlo based 
feasibility study. We generate B° —> d h° decays and process the events with 

detector simulation and reconstruction. Signal B candidates are selected. Signal and tagging 
B vertexes are reconstructed in order to obtain At, and the flavour of the tagging B meson 
is obtained. Finally, we perform an unbinned likelihood fit of the time-dependent Dalitz 
plot to obtain the value of 0i and its uncertainty. 


A. Monte Carlo Generation 


In order to test the feasibility of the method described above, we have developed an 

We first test 
In this case, 


3 . 


algorithm to generate Monte Carlo simulated data, based on EvtGen 
the generator by restricting the D —> decay to the KsP° channel. 














the formalism simplifies to the familiar Dcph 0 case, and the time-dependent decay rate 
(neglecting suppressed amplitudes), is given by 

g — I^I/ t bO 

V(At) = —--{1 + qS Dcph o sin(AmAf)} , (14) 

4t b o 

where the 6-flavour charge q is +1 (—1) when the tagging B meson is B° ( B °), and, within 
the Standard Model, Sn CP h° = ~VD C pVh°(— l) / sin(20i). For the CP odd decay D —> Kgp°, 
7 1d C p = ~ I; so f° r (KsP°)d 7r °’ ^D CP h° — ~~ sin(20!). Fig. QJshows generator level information 
for these decays. 



FIG. 3: At distributions for B° —> (Kstt + tt ) d h° (a) q = +1, (b) q = —1, (c) asymmetry. 


We next implement three body D decays into our generator. The amplitude of the 
D° —> K s tt + tt~ decay is described by a coherent sum of two-body decay amplitudes plus 
non-resonant part: 


f{m 


2 

Ks 


> m Kg«* 


N 

3 = 1 






Ks-k' 


' i m K s n+ 


) + be if> , 


(15) 


where N is the number of resonances, A, 


( m K s n+’ m Ks* 


+ h 


and o.j are the matrix element, 


amplitude and phase, respectively, for the j-th resonance, and b and [3 are the amplitude and 
phase for the non-resonant component. For further details, see H and references therein. 
Table □ describes the set of resonances we use in the decay model of our generator, which is 
similar to that in the CLEO measurement 0- Fig. El shows the Dalitz plot distribution for 
the D° —> KgTi + Ti decay generated according this model. 


For further confirmation of the operation of our generator, we look at the generator level 
time-dependent Dalitz plot. We generate B° —*• (K s tt + 7t~) D h° decays using 2 (f>i = 47°. In 
Fig. El we show the invariant mass distributions of the D decay daughters for events with 
q — — 1, and compare those for events with At greater than tb°/2 with those for events 
with At less than —tbo/2. Events with \At\ < tbo/2 or q — +1 are not shown. We see 
clear differences in the two invariant mass distributions; in particular we see more events 
with positive than negative At in the p° region of the 7r + 7r~ invariant mass distribution, as 
expected from Fig. El 

Since we are concerned with the feasibility of studying these modes at B factory exper¬ 
iments, we use the software of the Belle collaboration to perform simulation of the Belle 
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TABLE I: List of resonances used for D° —> I\stt + tt decay simulation. 


Resonance 

Amplitude 

Phase (°) 

/i*(892)+7r- 

1.418 

170 

/q(1430)+vr- 

1.818 

23 

/F 2 *(1430) + 7r- 

0.909 

194 

/t*(892)-7r+ (DCS) 

0.100 

341 

K S P° 

0.909 

20 

K s u 

0.034 

134 

A s ./o (980) 

0.309 

208 

A' 5 /o(1370) 

1.636 

105 

K s f 2 ( 1270) 

0.636 

328 

/i5'7r + 7r _ non-resonant 

1.0 

0 



FIG. 4: D° —> /i57r + 7r decay Dalitz plot. 



FIG. 5: Generator level invariant mass distributions of D decay daughters produced in the B° —> 
(Kstt + tt~) d h° decay chain. Events are generated with 2<f>i = 47°, and only events with q = —1 
(tagging B decays as B°) are shown. The dashed histograms show distributions for events with 
At > t b o/2, the solid histograms show those for events with At < —t b o/2. 
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detector and to reconstruct candidate events. The Belle detector is a large-solid-angle mag¬ 
netic spectrometer that consists of a silicon vertex detector (SVD), a central drift cham¬ 
ber (CDC), aerogel threshold Cerenkov counters (ACC), time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter (ECL) located inside a superconducting solenoid 
coil that provides a 1.5 T magnetic field. An iron flux-return located outside of the coil is 
instrumented to detect K Q L mesons and to identify muons (KLM). The detector is described 
in detail elsewhere [21]. The detector simulation is based on GEANT j|22]. Belle is in¬ 
stalled at the interaction point of the KERB asymmetric-energy e + e _ (3.5 GeV on 8 GeV) 
collider (23[. KEKB operates at the T(4S') resonance {\fs = 10.58 GeV) with a peak lumi¬ 
nosity that exceeds 1.5 x 10 34 cm~ 2 s _1 . The asymmetric energy allows At to be determined 
from the displacement between the signal and tagging B meson decay vertices. 


B. Event Reconstruction 


We reconstruct the decays B° —> (K S 7t + 7t~) d h° for h° = i r °,?7 and to. We use the 
subdecays Kg —> 7 r + 7 r _ , 7 T° —> 77 , 7 —> 77 , 7 r + 7 r _ 7 r° and u —* 7 r + 7 r _ 7 r°. The reconstruction, 
including suppression of the dominant background from e + e _ —> (Q = u > d, s, c ) continuum 
processes, is highly similar to that in related Belle analyses |24j. The properties of the 
background events are studied using generic BB and qq Monte Carlo. Our studies allow us 
to estimate the number of signal and background events to expect from a given data sample 
(we use the data sample of 253 fb -1 , containing 275 million BB pairs, collected with the 
Belle detector before summer 2004 as our baseline). The results are summarised in Table ITT! 


TABLE II: Detection efficiency, expected numbers of signal (Wig) and background (W,kg) events 
and signal purity for the B° —> (Ks'K + ir-) D h° final states. The expected numbers of events are 
based on the Belle data sample of 253 fb -1 . 


Process 

Efficiency (%) 

^sig 

-Vbkg 

Purity 

Dir 0 

8.1 

118 

49 

71% 

Duj 

3.9 

49 

8 

86% 

Dq 

4.3 

47 

15 

76% 

Sum 


214 

72 

75% 


For our further studies, we use only the Dn° mode, for which the expected number of 
signal events is the largest. In our pseudo-experiments, described below, we use numbers 
of signal and background events (300 and 100 respectively) which are rounded up from the 
totals in Table El as we expect some improvement is possible due to optimization of the 
selection for this analysis. 

The signal B meson decay vertex is reconstructed using the D trajectory and an in¬ 
teraction profile (IP) constraint. The tagging B vertex position is obtained with the IP 
constraint and with well reconstructed tracks that are not assigned to signal B candidate. 
The algorithm is described in detail elsewhere Si- 


Tracks that are not associated with the reconstructed B° —> (KsTT + TT~) D h° decay are 
used to identify the 6 -flavour of the accompanying B meson. The tagging algorithm is 
described in detail elsewhere [2fj]. We use two parameters, q and r, to represent the tagging 
information. The first, q , has the discrete value +1 (—1) when the tagging B meson is more 
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likely to be a B° ( B° ). The parameter r corresponds to an event-by-event flavour-tagging 
dilution that ranges from r = 0 for no flavour discrimination to r = 1 for an unambiguous 
flavour assignment. 

We divide the (pi = [0° : 180°] range into 18 points in steps of 10°. For each point 
we perform 30 pseudo-experiments with data samples consisting of 300 reconstructed Dn° 
events. We add 100 background events to each sample, where the background is modelled 
by B° —> D°h°, with uniform phase space decay D° —» Kstt +/ k~. 

For each pseudo-experiment, we perform a unbinned time-dependent Dalitz plot fit. The 
inverse logarithm of the unbinned likelihood function is minimized: 


- 21ogL = -2 


n 

E : 

i =1 


logp(m 2 , i ,m 2 _ i , Ati) — log / p(m 2 + , m 2 _, At)dm 2 + cim 2 _dAt 


D 


(16) 


where n is the number of events, m 2 +ll m 2 _ i and Ati are the measured invariant masses of 
the D daughters, and the time difference between signal and tagging B meson decays. The 
function p{;m 2 + , m 2 _, At) is the time-dependent Dalitz plot density, which is based on Eqs. 0 
and 0 including experimental effects such as mistagging and At resolution - we use the 
standard Belle algorithms to take these effects into account. The background component is 
also introduced into p. 




<l> 


1 input value 


<l> 


1 input value 


FIG. 6: (Left) average (pi fit result and (right) statistical error, as functions of the input value. 


Thus, for each input value of (pi we obtain fitted results from 30 pseudo-experiments. 
From the means and widths of the distributions of these results we obtain the average (j)\ fit 
results and estimates of their statistical errors. These results are shown in Fig. 0 We find 
the fit results are in good agreement with the input values, and the expected uncertainty on 
(pi is around 25°. 

To look for tails in the distributions, we also study larger ensembles of pseudo-experiments 
for two (pi input values: 23.5° and 66.5°, which correspond to sin(20i) = 0.73. We have 
performed this study both for the numbers of events described above, corresponding roughly 
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FIG. 7: Fit results for 0i. The solid (hatched) histograms correspond to the input value 20! = 47° 
(20i = 133°). The left (right) plot corresponds to a data sample roughly equivalent to 250 ftT 1 
(500 fb' 1 ). 


to 250 fb _1 (for which we perform 250 pseudo-experiments for each input value of 0x), and 
for numbers twice larger (hence 500 fb _1 , for which we perform 125 pseudo-experiments). 
Fig-IUshows the distribution of the fit results. We do not observe any pathological behaviour, 
demonstrating that this method can indeed be used to distinguish the two solutions for 
sin(20i) with sufficiently large data samples. 

We have tested for possible bias in the method due to neglect of the suppressed amplitudes 
(Eqs. C3 and 0. Due to the smallness of r Dn o compared to the B° — B° mixing effect, we 
expect any such bias to be small, and indeed we find it to be smaller than 1%. 

As noted above, this method is highly similar to that used to extract 03, using B ± —*■ 
DK ± followed by multibody D decay BUS- A significant complication arises in that case 
due to uncertainty in the D decay model, and we expect this will also affect the B Q —> Dh° 
analysis. However, the time-dependent analysis does not suffer due to the smallness of the 
ratio of amplitudes, and therefore we expect that the model uncertainty may be smaller. 
Furthermore, a number of methods have been proposed to address the model uncertainty 
(for example, using information from CP tagged D mesons which can be studied at a cr 
factory, such as CLEO-c), and this analysis can also take advantage of any progress in that 
area. 


IV. CONCLUSION 

We have presented a new method to measure the Unitarity Triangle angle 0i using am¬ 
plitude analysis of the multibody decay of the neutral D meson produced in the processes 
B° —>■ Dh°. The method is directly sensitive to the value of 20! and can thus be used to 

resolve the discrete ambiguity 20! <-> 7r — 20!. The expected precision of this method has 

been studied using Monte Carlo simulation. We expect the uncertainty on 0! to be about 
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25° for an analysis using a data sample of 253 fb 1 
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